High Order Embedded Runge-kutta Scheme for Adaptive Step-size Control in the Interaction Picture Method
نویسندگان
چکیده
The Interaction Picture (IP) method is a valuable alternative to Split-step methods for solving certain types of partial differential equations such as the nonlinear Schrödinger equation or the Gross-Pitaevskii equation. Although very similar to the Symmetric Split-step (SS) method in its inner computational structure, the IP method results from a change of unknown and therefore do not involve approximation such as the one resulting from the use of a splitting formula. In its standard form the IP method such as the SS method is used in conjunction with the classical 4th order Runge-Kutta (RK) scheme. However it appears to be relevant to look for RK scheme of higher order so as to improve the accuracy of the IP method. In this paper we investigate 5th order Embedded Runge-Kutta schemes suited to be used in conjunction with the IP method and designed to deliver a local error estimation for adaptive step size control.
منابع مشابه
Embedded Runge-Kutta scheme for step-size control in the interaction picture method
When solving certain evolution type PDE such as the Schrödinger equation, the Interaction Picture method is a valuable alternative to Split-Step methods. The Interaction Picture method has good computational features when used together with the standard 4th order Runge-Kutta scheme (giving rise to the RK4-IP method). In this paper we present an embedded Runge-Kutta scheme with orders 3 and 4wit...
متن کاملSPIP: A computer program implementing the Interaction Picture method for simulation of light-wave propagation in optical fibre
The computer program spip is aimed at solving the Generalized Non-Linear Schrödinger equation (GNLSE), involved in optics e.g. in the modelling of light-wave propagation in an optical fibre, by the Interaction Picture method, a new efficient alternative method to the Symmetric Split-Step method. In the spip program a dedicated costless adaptive step-size control based on the use of a 4th order ...
متن کاملComparison of adaptive step-size control strategies for solving the Generalised Non-Linear Schrodinger Equation in optics by the Interaction Picture method
This report is devoted to the study and to the comparison of methods for estimation of the local error and for adaptive step-size control when solving the generalised nonlinear Schrodinger equation (GNLSE) in optics by the Interaction Picture (IP) method or by the Symmetric Split-Step (S3F) method. Namely, we propose and study the use of an embedded Runge-Kutta method to solve the nonlinear pro...
متن کاملEmbedded implicit Runge-Kutta Nyström method for solving second-order differential equations
An embedded diagonally implicit Range-Kutta Nystrom (RKN) method is constructed for the integration of initial value problems for second order ordinary differential equations possessing oscillatory solutions. This embedded method is derived using a three stage diagonally implicit Runge-Kutta Nystrom method of order four within which a third order three stage diagonally implicit Runge-Kutta Nyst...
متن کاملEmbedded Diagonally Implicit Runge - Kutta Algorithms on Parallel Computers
This paper investigates diagonally implicit Runge-Kutta methods in which the implicit relations can be solved in parallel and are singly diagonalimplicit on each processor. The algorithms are based on diagonally implicit iteration of fully implicit Runge-Kutta methods of high order. The iteration scheme is chosen in such a way that the resulting algorithm is ^(a)-stable or Z,(a)-stable with a e...
متن کامل